Motional dressed states in a Bose condensate: Superfluidity at supersonic speed 
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We present an exact analytic solution of a nonlinear 
Schrodinger field interacting with a moving potential (obsta- 
cle) at supersonic speed. We show that the field forms a stable 
shape-invariant structure localized around the obstacle — a 
dressing effect that protects the field against excitations by 
the obstacle. 

PACS numbers: 03.75.Fi, 67.57.De 

One of the most intriguing phenomena in superfiuid is 
frictionless flow below a critical velocity v c Q . This phe- 
nomenon, first observed in liquid helium, demonstrates 
key features of collective behavior in macroscopic quan- 
tum coherent systems. Recently, experimental evidence 
of the existence of a critical velocity in a Bose-Einstein 
condensed gas has been found Q . It is known that atoms 
in such a dilute system interact weakly with each other, 
and so it justifies the use of the nonlinear Schrodinger 
equation (NLSE) to study the quantum dynamics of the 
condensate. Indeed, direct numerical solutions of the 
NLSE have indicated a distinct transition from superfiuid 
flow to resistive flow at a critical velocity ^Q. Recent 
studies have also addressed dynamic features in solutions 
of the NLSE, such as phonon and vortex emissions 
in order to understand the dissipation mechanisms. 

As an obstacle moves through a Bose condensate, the 
regime of supersonic speed is often considered as a dis- 
sipation domain where frictionless motion disappears. 
This is understood from Landau's argument that energy- 
momentum conservation forbids phonon excitations un- 
less the obstacle's velocity is at least the speed of sound 
c, i.e., the critical velocity v c equals c for systems with a 
pure phonon spectrum In fact, for realistic systems 
obeying the NLSE, v c is less than c because of the re- 
duced fluid densities appearing at the boundary of the 
obstacle ]^-|5|,|To| . Therefore the speed of sound seems 
to be an upper limit beyond which dissipation occurs 
generally. However, we shall show that there exist in- 
teresting exceptions at least in one dimensional systems. 
Such exceptions exhibit rich dynamic features resulted 
from coherent field-obstacle interactions, which cannot 
be captured by the concept of critical velocities alone. 

In this paper we present an exact analytic solution of 
a nonlinear Schrodinger field interacting with a moving 
repulsive potential (obstacle) in one dimension. Our so- 
lution represents a family of motional dressed states in 
which the field organizes itself as a stable and shape in- 
variant structure localized at the moving obstacle. Once 



the dressing is fully developed, there is no energy transfer 
from the obstacle to the field. We show that this hap- 
pens when the velocity of the potential is greater than the 
speed of sound in the field. Therefore dressed states for- 
mation is a novel coherent feature that can be maintained 
in a condensate at high speeds. The existence of dressed 
states at supersonic speed requires a specific matching 
between the field and the potential. Our solution pro- 
vides a prescription for the matching potentials. This 
opens possibilities of coherent control of local properties 
of a Bose condensate. 

We begin by writing down a nonlinear Schrodinger 
equation for a field ^(x, t) under the influence of a mov- 
ing potential U{x — vt), which is well localized and re- 
pulsive. As usual, the NLSE in a dimensionless form is 
given by 

ijU(x,t) = -i j^OM) + Wx,t)\ a 9(x,t) 

+U(x-vt)^(x,t) -V(x,t). (1) 

This equation can be used to determine the evolution 
of the condensate wavefunction in one dimension at zero 
temperature. If the condensate has a chemical potential 
fi and a particle mass m, then the position x and the 
time t are in units of %/y/rful and h/n, respectively. In 
Eq. (I), the nonlinear term is positive for repulsive self 
interactions, and the speed v is taken to be positive for 
definitcness. We assume that the field is uniform (i.e., 
at rest) as x tends to ±oo. This imposes the boundary 
condition 

*(x -► ±oo,*) = e i7± (2) 

at any finite time *. Here j± are constants which account 
for possibilities of having different phases at the two in- 
finities. The choice of unit modulus is a convenient scale 
which gives the speed of sound c = 1 for undisturbed 
systems. 

We define a dressed state of the system by a solution of 
(1) in the shape invariant form: ^(x, *) = ip[x — vt) sub- 
jected to the boundary condition (2). In other words 
dressed states are stationary states in the co-moving 
frame of the potential. For v in the subsonic domain 
v < 1, Hakim has reported stationary states which are 
of the form of a dip surrounding the obstacle pp[ . The 
dip appearance is understood because of the repulsive 
character of the potential. However, for v greater than 
the speed of sound, no known steady states exist. We 
now describe a method to determine dressed states in 
this regime. 
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By assuming the solution in the form: iP(t) 
r(r)e' e ' T - ) , where r = x — vt, Eq. (1) requires 



- rr = [r A -r)+ U{r)r - vrO 



- (rij + 2f 



(3) 
(4) 



Here the derivatives are taken with respect to r, i.e., 
r = dr/dr,r = d 2 r/dr 2 , etc. If we consider r as a kind 
of time, then equations (3) and (4) are equivalent to the 
motion of a classical particle interacting with a central 
force and a time-dependent magnetic field described by 
U(t). Therefore a dressed state corresponds to a trajec- 
tory matching the boundary conditions (2). This requires 
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±00) = 1, 
±00) = 7-1- . 



(5) 
(6) 

- r 2 9, we see that 
K must equal v for desired trajectories. In this way we 



From the constant of motion K = 



have 8 — v(l 



and so Eqs. (3) and (4) are reduced 



to a one-dimensional trajectory problem 
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U(r)r. 



(7) 



Once the function U (r) is given, we may search for a 
solution by finding out initial conditions r(0) and f (0) 
that make r(r) approach one asymptotically. Although 
the search can be performed systematically via numerical 
means, the existence of a solution is not a guarantee for 
arbitrary functions U and speeds v. 

We obtain dressed state solutions analytically by 
adopting an inverse approach. First we construct an r(r) 
that satisfies the boundary condition (5), then we deter- 
mine the corresponding U(t) from Eq. (7). The require- 
ment is that the U (r) must be well localized and every- 
where repulsive, i.e., [/'(£> 0) < and U'(£ < 0) > 0. 
This approach tells us how one should design the shape 
of the U in order to have a specified form of 'J. The 
matching of U(t) and r(r) yields the exact dressed state 
solutions. We note that in experimental situations [Q, U 
can in fact be an adjustable function controlled by the 
intensity of a detuned laser. 

We give a family of exact solutions in the form of a 
solitary wave: 



r(r) = [l- -sech 2 (/3r) 

L v 

Q( T ) = ~ tanh /3r. 

P 



-1/2 



(8) 
(9) 



The two positive parameters a and (3 characterize the 
height and inverse width of the field \& near the moving 
potential. The exact form of U(t) that matches Eq. (8) 
is given by, 

tt, \ asech 2 (/3r) , ,i, ns 

U{t) = K — 2 l A i + A 2 sech 2 (/3r) 

2 [w-asech 2 (/3r)] 

+ A 3 sech 4 (/3r) + A 4 sech 6 (/37-)} (10) 



where Ai = 2v (v 2 + 1 - l) , A 2 = 2a + a/3 2 - 3f3 2 v 
5av 2 , A3 = 4a 2 v and A4 = - 



We find that if v satisfies 



(v - a) 3 



1 > 2/3 2 , 



(11) 



then U(t) in Eq. (10) is repulsive as required. Indeed, 
v has to be supersonic in order to satisfy the inequality 
(11). Therefore we have provided an explicit solution of 
dressed states at supersonic speeds. 

In Fig. 1, we illustrate the general shapes of the match- 
ing potentials U(t) as a function of parameters a and [3. 
We plot a characteristic curve defined by the equality sign 
in (11). Therefore the region under the curve corresponds 
to purely repulsive potentials. These potentials have a 
maximum barrier height at r = 0. Above the curve, the 
second derivative U"(0) is positive which corresponds to 
an attractive force on particles near the center of U. We 
shall not discuss the solutions in such a domain, since the 
potentials are not purely repulsive. We remark that at 
a higher speed v, the characteristic curve takes a higher 
position on the a — (3 plane, allowing a wider range of 
parameters for the solution (8-9). We also note that 
U(t) takes a simpler form, U{t) w wasech 2 (/3r) , when 
v > a,fl. 

There are interesting features in our solution: (I) A 
dressed state for a repulsive U at v > 1 must come with 
a positive a. This means that the field concentrates at 
the center of U, even though U itself is repulsive. Such 
a 'counter-intuitive' behavior does not occur for dressed 
states in the subsonic speed domain. It is worth noting 
that the family of solutions Eqs. (8-10) can be general- 
ized for subsonic speeds. In that case, a is negative which 
recovers the expected dip at the center of the potential. 
(II) The height of the matching potential U is bounded 
by 



2*7(0) < (w 2 -3u 2/3 + 2). 



(12) 



This inequality sets an upper limit of U (0) for the exis- 
tence of the solution (8-9). From the right side of Eq. 
(7), we see that a large U term would lead to an un- 
bounded r(r) unless the v 2 term can make a balance. 
Hence for speeds v which are close to (but still greater 
than) 1, the matching potential must be weaker accord- 
ingly. It is also apparent that there are no dressed states 
for impenetrable potentials. 

In order to investigate the dynamical consequences of 
dressed states, we solve numerically the time-dependent 
NLSE (1) for ^(x,t) which starts from the ground state 
\f(a;,0) = 1. The potential U(x — vt) is turned on sud- 
denly at t = 0. In Fig. 2a, we plot the field energy, 

AE(t) = \jdx (|V*(a;,i)| 2 + |*(x,t)| 4 ) -£7(0) rela- 
tive to the initial field energy £7(0) as a function of time. 
This quantity tells us how much field energy is gained 
from the moving potential if it is switched off at a time t. 
In Fig. 2a we see that there is a positive gain of the field 
energy at early times. After a characteristic time of or- 
der 2/3 -1 , the energy eventually becomes a constant. In 



2 



this steady state regime, the moving potential does not 
'feel' any the reaction force from the field even though U 
is repulsive. 

The emergence of a steady state indicates the forma- 
tion of a dressed state. This is explicitly shown in Fig. 
2b where are plotted as a function of x at different 
times. The arrows in the figure indicate the positions 
of the moving potential. We see that the field develops 
a shape invariant wavepacket form at the potential, i.e., 
the potential is dressed. By comparing both the phase 
and the amplitude of the wavepacket with the analytic 
expressions Eqs. (8) and (9), we found a very good agree- 
ment between them. At regions far away from the poten- 
tial, \& exhibits features created by the sudden motion at 
t = 0. We see that the left-most wavefront propagates 
at the speed of sound, and there is an oscillatory pattern 
developing in a depression region near the middle. These 
features are partly responsible for the transient energy 
gain shown in Fig. 2a. A closer look at the depression re- 
gion at longer times (not shown) indicates that the short 
wavelength structure becomes grey solitons traveling at 
subsonic speeds pT[ . 

We emphasize that the usual picture of critical ve- 
locities does not apply to dressed states. The velocity- 
dependent requirements here are specified by the inequal- 
ities (11) and (12), which do not impose an upper limit 
for v. Instead, v has to be sufficiently high in order 
to achieve dressing for a given U JT^j. We have calcu- 
lated the energy gain using the same potential shown in 
Fig. 2a but at different velocities. For velocities near 
the speed of sound such that (12) is strongly violated, 
we see no dressed state formation. In those cases the 
field energy keeps growing via repeated emissions of grey 
solitons from the potential as observed in Ref . [1(J . 

Finally, we want to point out the stability of our so- 
lution. Since the initial condition used in Fig. 2b is far 
away from the dressed state solution, the convergence of 
the field precisely into the analytical form (8) and (9) is 
an evidence of dynamic stability of dressed states. We 
have performed direct stability tests based on collective 
excitation frequencies in the co-moving frame. Our nu- 
merical results confirm the dynamic stability of the solu- 
tions Jl|. The issue of the sensitivity to the shape of U 
has also been examined. By replacing U with a gaussian, 
we repeat the numerical calculations as in Fig. 2b. As 
long as the gaussian's width and height are well within 
the domain described by inequalities (11) and (12), we 
see the emergence of dressed states. Therefore dressed 
state formation is not sensitive to the details of U. 

To conclude, we have presented a family of motional 
dressed states as exact analytic solutions of a nonlin- 
ear Schrodinger field interacting with a moving potential. 
The formation of a dressed state enforces frictionless mo- 
tion, and so it preserves a key feature in superfluid in 
the supersonic regime. We stress that our method is 
not limited to the solution family (8-9) and potentials 
in purely repulsive forms. With our approach, one can 
construct wide varieties of solutions of the field and the 



corresponding matching potentials, as long as the sta- 
bility and physical constraint for the potential are al- 
lowed. This opens interesting possibilities of controlling- 
local densities of a condensate dynamically. Our study 
is for one-dimensional systems; we speculate that similar 
dressing features appear in Bosc condensates in quasi- 
1D configurations, such as condensates trapped in a long 
cigar shape or in a ring, within characteristic time scales. 
It remains an important open question whether dressing 
effects appear in two or higher dimensional systems. We 
hope to address this issue in the future. 
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M -{ -r 2 -D*-2\i>\ 2 ) 

where D = —\d 2 , + ivd x /. The system is dynamically 
stable if all eigenvalues of M are real. 
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FIG. 1. Parameter space of the solutions at a supersonic 
speed v = 2, with positive a and 13. The characteristic curve 
is defined by the equality sign of (11). The region under the 
curve is associated with purely repulsive potentials U . The 
shapes of U at the two dots are representative examples above 
and below the curve 
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FIG. 2. (a) The increase of energy as a function of time 
with the parameters: v = 2, a = 0.2 and j3 = 0.05. The 
potential is shown in the inset, (b) Time evolution of the 
modulus of ^ at different times. The curves are shifted up 
by 0.1 successively in order to display the evolution. The 
arrows indicate the position of the center the potential at the 
corresponding time. 
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